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Abstract 

The quantum mechanical two-body problem with a central interaction on the 
sphere S" is considered. Using recent results in representation theory an ordinary 
differential equation for some energy levels is found. For several interactive potentials 
these energy levels are calculated in explicit form. 

PACS numbers: 03.65.Fd, 02.40. Vh, 02.40. Ky. 

Mathematical Subject Classification: 43A85, 22E70, 57S25, 70G65. 

1 Introduction 

The history of mechanics on constant curvature spaces started more than one and a half 
century ago. 

The analog of Newton (or Coulomb) force for the hyperbolic space H'^ was proposed 
already by founders of the hyperbolic geometry N.I. Lobachevski (in 1835-38) fl] and 
J. Bolyai (between 1848 and 1851) 2 as the value F(p), which is inverse to the area of 
the sphere in H'^ of radius p with an attractive body in the center. 

The analytical expression for the Newtonian potential in the space H'^ was written 
in 1870 by E. Schering 21 (see also his paper of 1873), without any motivation and 
references to N.L Lobachevski and J. Bolyai. 

In 1873 R. Lipschitz considered a one-body motion in a central potential on the sphere 

IS]. He knew that the central potential Vc satisfies the Laplace equation on S'^. However 
due to some reason he preferred to consider another central potential V{p) ^ sin^^(p/i?), 
where p is a distance from the center and i? is a curvature radius. He calculated the 
general solution of this problem through elliptic functions. 

In 1885 W. Killing found the generalization of all three Kepler laws for the sphere 
[H]. He considered the attractive force as an inverse area of a 2-dimensional sphere in 
as N.I. Lobachevski and J. Bolyai did before. In the next year these results was published 
also by C. Neumann in Their expansion onto the hyperbolic case was carried out in 
H. Liebman paper 8 in 1902 and later in 1905 in his book on noneuclidean geometry |^. 
Note that he started from ellipses in or and derived a potential in such a way that 
the first Kepler law would be valid. He derived also the generalization of the oscillator 
potential for these spaces from the requirement that a particle motion occurs along an 
ellipse with its center coinciding with the center of the potential. 

Also in the same paper [HI W. Killing proved the variable separation in the two-centre 
Kepler problem on the sphere S", which implies the integrability of this problem. 

The well-known Bertrand theorem ^10 states that up to an arbitrary factor there are 
only two central potentials in Euclidean space that make all bounded trajectories of a 
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one-body problem closed. In spaces S^, also there are only two potentials Vc and Vo 
with this property. It was proved by H. Liebman in 1903 ^Jj, see also P|. 

One can consider the classical mechanics in spaces of constant curvature as a prede- 
cessor of special and general relativity. After the rise of these theories the above papers 
of E. Schering, W. Killing, H. Liebmann were almost completely forgotten. Note that the 
description of a particle motion in central potentials in spaces and was shorten in 
the second and the third editions of the H. Liebman book 9 w.r.t. the first edition in 
favor of special relativity. 

Similar models attracted attention later from the point of view of quantum mechanics 
and the theory of integrable dynamical systems. This leads to the rediscovery of results 
described above in many papers. Note however that almost forgotten results of E. Schering, 
W. Killing and H. Liebmann were described in the survey |12|. 

Quantum mechanical spectral problem on the sphere S'^ for potential Vc (Coulomb 
problem) was solved by E. Schrodinger in 1940 by the factorization (ladder) method, 
invented by himself A.F. Stevenson in 1941 solved the same problem using more 

traditional analysis of the hypergeometric differential equation ^3] (see also L. Infeld 
result in 1941 JHl)- L. Infeld and A. Schild in 1945 solved a similar problem in the space 
H3 [m (see also [HI). 

The connection of the Runge-Lenz operator for the quantum Kepler problem in 
with the Schrodinger ladder method was discussed by A.O. Barut and R. Wilson in |18j . 
In A.O. Barut, A. Inomata and G. Junker solved the Kepler problem in and 
using the functional integration. 

In papers [20], V.S. Otchik considered the one particle quantum two center Coulomb 
problem in S'^ and found a coordinate system admitting the variable separation. The cor- 
responding ordinary differential equations are Heun's ones. 

In [25] - there was developed an algebraic approach to one particle problems for 
potentials Vc and Vo in spaces S",H". 

Transformations between the Coulomb-Kepler and oscillator problems existing in the 
Euclidean space were generalized for the sphere in j2.7J . 

In 28 A. A. Bogush, Yu.A. Kurochkin and V.S. Otchik considered the Coulomb scat- 
tering in the space H^. 

The two-body problem with a central interaction in constant curvature spaces S" and 
H" considerably differs from its Euclidean analog. The variable separation for the latter 
problem is trivial, while for the former one no central potentials arc known that admit a 
variable separation. 

The classical two-body problem with a central interaction in constant curvature spaces 
was considered for the first time in |29j . Its Hamiltonian reduction to the system with two 
degrees of freedom was carried out by explicit coordinate calculations. For some potentials 
there was proved the solvability of the reduced problem for an infinite period of time. 

In (SOI there was studied the self-adjointness of the quantum two-body Hamiltonian in 
spaces and and were found in explicit form some its infinite energy level series for 
the sphere S^, corresponding to some central potentials. 

Simply connected constant curvature spheres S" and hyperbolic spaces H" are rep- 
resentatives of the class of two-point homogeneous Riemannian spaces (TPHRS). Such 
spaces are characterized by the property that any pair of points can be transformed by 
means of an appropriate isometry to any other pair of points with the same distance 
between them. Equivalently, these spaces are characterized by the property that the nat- 
ural action of the isometry group on the unit sphere bundle over them is transitive. The 
classification of TPHRS can be found in |31| . 

For a smooth manifold M endowed with a left action of a Lie group G denote by 
Diff(M) = DiSoiM) the algebra of G-invariant differential operators on AI with smooth 
coefficients. For a Riemannian manifold M let Ms be the unit sphere bundle over M. 
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Let Q be an arbitrary TPHRS, endowed with the action of the identity component of the 
isometry group for Q. 

In 1221 there was found a polynomial expression for the quantum two-body Hamiltonian 
H on Q through a radial differential operator and generators of the algebra DiS{Qs)- 
Coefficients of this polynomial depend only on the distance between particles. 

Algebras Diff(Q5) are noncommutative. A full set of their generators and correspond- 
ing relations^ was found in [55) . 

Let 21 be a set of Diff (Qg) generators presented in the expression for the Hamiltonian 
H. An every common eigenfunction of operators from 21 generates a separate spectral 
ordinary differential equation for the two-body quantum mechanical problem on TPHRS. 
The search for such common eigenfunction is not an easy problem. In low dimensions for 
Q — S^, Q = this problem was solved in [30] and [23] using an explicit description of 
S0(3) and S0(4) irreducible representations. The present paper deals with this problem 
for the general spherical case Q — S". The progress is reached using the results in 
representations theory of the algebras so(n,C) in [23 and pTH] 

The paper is organized as follows. Sections [3H are of a preparatory character. Section 
[21 and [31 contain basic facts on invariant differential operators on homogeneous spaces 
and regular representations of compact Lie groups respectively. In section ^ there is a 
description of the quantum two-body Hamiltonian on the sphere S" through a radial 
differential operator and generators Di, z = 0, 1, 2, 3 of the algebra Diff (S^). 

Sections [HI and [Hi form a main part of the paper. In section [3 we calculate actions 
of operators Di, i = 0,1,2,3 in a corresponding functional space and find all common 
eigenvectors ipD for operators Dq,Di,D2 and optionally D3. Using these eigenvectors 
we derive in section [Sj a separate ordinary differential equation of the second order for 
a radial part of a two-body eigenfunctions. For Coulomb and oscillator potentials this 
differential equation is Fuchsian and we consider its reducibility to the hypergeometric 
one using the rational change of an independent variable. This reduction is possible for 
some eigenvectors ip d that leads to an explicit form of some infinite energy level series for 
the two-body problem with Coulomb and oscillator potentials. 

A necessary information concerning complex orthogonal Lie algebras, self-adjoint Schro- 
dinger operators on Riemannian spaces and Fuchsian differential equations is collected in 
appendices IXHCl 

ACKNOWLEDGEMENT. The author is grateful to A.I. Molev for pointing out 
his papers ^Sj and 

2 Invariant differential operators on homogeneous 
spaces 

Here we shall briefly describe the construction of invariant differential operators on homo- 
geneous spaces [23- 

Let G be a Lie group of dimension iV and K be its subgroup of dimension N ~ £. 
Denote the corresponding Lie algebras as g and B. Suppose that the algebra g admits the 
reductive expansion 

= p®t, (1) 

for a subspace peg, i.e. Adi^ p C p. For a compact Lie group G such subspace p can 
always be constructed using the invariant integration on G. Let {ej)jLi be a base in g 
such that {ejYj^i is a base in p. 

^One relation for the quaternion projective space and its hyperbolic analog was calculated only in 
leading terms. 
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Let S'(p) be a symmetric algebra for the linear space p. The Ad/f-action on p is 
naturally extended to the Ad/f -action on S{p). The main result of the general theory [57] 
is that G-invariant differential operators on G/K are in one to one correspondence with 
the set S{p)^ of all Adif-invariant elements in S{p). 

Let I : p ^ S'(p) be an inclusion, U{q) be the universal enveloping algebra for g and 
A : ^(p) — > U{q) be a linear symmetrization map, defined on monomials by the formula 

■ ■ ■ ■ ■ = ^ • • ■ ■ • ev(fc) :i i^ij ^ i ^k, 

where e*^. := i{ei-) and 6^ is the full permutations group of k elements. Evidently 

A : 5(p)^ ^ C/(0)^, 

where [/(fl)''^ is the set of all Adi^-invariant elements in J7(g)^. 

Let P{ei, . . . jCn) be a polynomial depending on noncommutative elements. Denote 
by Ci the left invariant vector field on G, corresponding to the element a € g = Tf,G : 

5exp(tei), 5 e G. 

t=o 

Then Z?p := P{ei, . . . , ejv) is a left invariant differential operator on G. 

Functions on the homogeneous space G/K are in one to one correspondence with 
functions on the group G that are invariant w.r.t. right iiT-shifts. For P(ei, . . . , ca?) € 
U{g)^ the differential operator Dp, acting on such functions, can be considered as a G- 
invariant differential operator on the space G/K and every such operator can be uniquely 
represented in the form 

(A(Po))(gi,...,e£)), 

for some Pq G 5'(p)^. 




3 Regular representations of compact Lie groups 

Let G be a compact connected Lie group and /Lt be a biinvariant positive measure on 
G, unique up to an arbitrary factor jHEl- Let £^(G, /x) be a Hilbert space of measurable 
complex valued functions on G, square integrable w.r.t. the measure fj,. Define two unitary 
left representations of G in the space C^{G, /i). The left regular representation T' acts by 
the left shifts 

(r^/) (5) = fiq-'g), q,geGje ii\G, 

and the right regular representation T'' acts by the right shifts 

i^li) (5) = /(ot), g,5 e G,/ e ^^(G,/.). 

Evidently these representations are equivalent with the intertwining operator /(g) — > 
f{g~^)- It is well known that these representations expand into direct sums of finite 
dimensional unitary irreducible representations (irreps) . Each of these irreps is contained 
in or with a multiplicity equal to its dimension and an every linear irreducible 
representation of G is equivalent to an irreps from this sum [35], |4l)j. 

Let Ti be a full system of unitary irreps for G in spaces Ug, £ = 1,2,.... Choose in every 
Ui an orthonormal base (e^,/c)feLii d,i := dime Ui. Define matrix elements t^ ^ of operators 

by the equation T^eg^k t\_k{q)(^t,i or equivalently by t\j^{q) := {ef,^i,T^eg^k)un q ^ G. 
Since 

t\,ki9q)et,t = T^T^ee^k = ti,i9)tlkiqye,j, g,q^G 
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one has 

tiki9Cl)=ti,{9)ti,{q). (2) 

Therefore the subspace TZi^i C C^{G, ji), spanned by functions {t\ ■ {g))'^ti i is invariant un- 
der operators and the representation T^\-ji^ is equivalent to Ti. On the other hand, the 

formula (0) impHes that the subspace Cij C C'^{G,iJ,), spanned by functions (i^ 

is invariant under operators and the representation is again equivalent to Tg. 

The functions {t\ j{g))'l'^j^n I = 1,2, .. . form an orthogonal base in the space C?(G,ii) 
EHI, EHI, mi and 

||,» ||2 _ ^^jG) 

Thus the space 

:=07^,,, = 0A,, 
^=l J=l 

is invariant under representations T'' and T'. The representation T'" intermixes spaces 
Hi^j of representations T' and vise versa the representation T' intermixes spaces Tl-i^i of 
representations T''. The space C?{G,\x) of representations T'' and T' expands into irreps 
as follows 

c? (G, = r, = 7^,.. = 00 A, . 

i e i=i i j=i 

For a Lie subgroup K of the group G the subspace (G, if, /i) C C?[G, ^), consisting of 
functions invariant w.r.t. all right -ftT-shifts on G, is invariant w.r.t left G-shifts. Therefore 
there are only two possibilities: 

n £2 (G, K, fi) = Cij and £f j n (g, if, ^) = 0. 

The consideration above implies the following proposition. 

Proposition 1. Let 

ft := Ti n £2 (G, K, ^) , := 7^^,, n (G, fi) , := dime ^f.z- 



Evidently, the value di does not depend on i — l,...,dg. The representation ^ |^ is 

expanded into the direct sum of equivalent irreps in spaces Cf^, k — 1, . . . ,dg, which are 
among of Cg^j . On the other hand 

dt ^ 

= 07^£,„ 
1=1 

where the spaces TZi^i, i — I, . . . , dg are isomorphic to each other. 

4 Two-body Hamiltonian on the sphere S"^ 

Let S", n 5^ 2 be the n-dimcnsional sphere, endowed with the standard metric g of a 
constant sectional curvature R^^, R > and 
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be the corresponding Laplace-Beltrami operator, expressed through local coordmates, 
where 7 :— det \\gij\\. We start from the description of the two-body quantum Hamil- 
tonian on S" found in and gTj . 

The configurations space for the two-body system on S" is 

S" X S". (3) 

The Hamiltonian for this system is 

Hv = Ho + V=--^ A,^-^ A2+V{p), (4) 
2toi 2to2 

where Ai, i = 1,2 is the Laplace-Beltrami operator on the ith factor of Q and p be 
the distance between particles. It should be defined on some subspace Dom(ff ) dense in 
£^ (S" X S",x X x) to be a self-adjoint operator, where x is the measure on S" induced 
by the metric. In local coordinates x has the form: x = y/jdxi A ... A dxn- Note that the 
free Hamiltonian Hq is the Laplace-Beltrami operator for the metric 

52 := TOi7r*g -|- m2TT2g (5) 

on multiplied by —1/2, where 7r*(7 is the puUback of the metric g with respect to the 
projection on the i-th factor. 

Let G = SO(n -I- 1) be the identity component of the isometry group for the sphere S". 
One can consider SO(n -I- 1) in the standard way as 

SO(n + 1) = {Ae GL(n + 1,M)| AA^ = E, dctA = l) , 

where E is the matrix unit. The configuration space © is endowed with the diagonal 
G-action and the differential operator 0) is G- invariant. 

Let K = SO(n — 1) be a subgroup in SO{n + 1) with elements of the form 

A e SO(n- 1). 

Up to a manifold of dimension n, consisting of antipodal points, the configuration space 
© can be represented as the direct product 

/ X {G/K) , (6) 

where / = (0, ttR) and the factor space G/K is G-homogeneous w.r.t. left shifts [32]. The 
space G/K is isomorphic to the unit sphere bundle over S" . 

The Lie algebra g = so(n-l-l), corresponding to the group G, consists of skew-symmetric 
matrices. Let E^j be the matrix of the size (n -1- 1) x {n + 1) with the unique nonzero 
element equals 1, locating at the intersection of the fc-th row and the j-th column. Choose 
the base for the algebra g as: 

= - Ejk, 1 s; fc < j n + 1. 

The algebra g admits the reductive expansion where the subspace p is spanned by 
elements 

^-ife, 2 s; fc n + 1, *2fc, 3 fc s: n + 1. 

In the general case n ^ 4 generators of the commutative algebra S{p)^ can be chosen 
as: 

n-\-l n+1 n+1 

fz — 3 A, — 3 — 3 



/ E2 










\ 




, 
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In the case n — 3 there is the additional generator 
In the case n — 2 the group K is trivial and generators of S{p)^ — S{p) are simply 

*t2, *l3> *23- 

In all cases we shall consider elements 

n+l n+1 ^ n+1 

Do = -vfia, A - = 51 *2fc, ^3 - -2 E {^i'^' *2fe} , (7) 

from U{q) as invariant differential operators on the space G/K, where {•,•} means an 
anticommutator. The commutative relations for differential operators lO are (see 

[Do, Di] - -2D^, [Do.D^] - 2i?3, [Do, D^] = Di - D2, [Di,D2] - ~2{Do, D3}, (8) 
[D,,D,] = -{Do, D,} + i-^)i--^) n„ [D2,D,] = {Do, D,} - ili^lKli^^,. 
For n — S the additional operator 

□ := ^ ({*13> *24} - {^-w, ^'23}) 

lies in the centre of the algebra DiSG{G/ K). 

Define a new coordinate r on the interval / by the equation 

r = tan(^) , r eM+ (0,oo). 

Results from and imply the following theorem. 

Theorem 1. The quantum two-body Hamiltonian on the sphere S" can be considered as 
the differential operator 

SmR'^r''-^ dr° \ + r2)"-2 dr J 2mim2R^ " ^ ^ 

4mim2it-^r" ^ [_ or (l + r^)" 2 
on i/ie space M_|_ x G, where 

mim2 

m :— , (10) 

mi + m2 

a parameter a G (0, 1) is arbitrary, P :— I — a, and 
(1 + r2)2 

^ = no (mi cos^(2Q;arctanr) + 7712 cos^(2/?arctanr)) , 

4miTO2^ 
(1 + ^2)2 

B = TT-ir (mi sin(4a arctanr) — m2 sin(4/3 arctanr)) , 

8mirn2^ 
n j,2\2 

C = 7^ (nil sin2(2Q; arctanr) + m2 sin^ (2/3 arctanr)) . 

^mim2R r'^ ^ ' 

The domain for operator ^ is dense in the space (M+ x G,K,rj), consisting of all 
complex valued square integrahle K -invariant functions on x G, with respect to right 
K-shifts and the measure 

r"-''-dr 

dv = -r, ^TT' ® an = dv (i) da, 

(1 + r'^)" 

where n is a biinvariant measure on G, unique up to a constant factor. 
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In the following we choose the parameter a in such a way that mia — m2/3 — 0, i.e. 

m2 „ mi 



mi + m2 mi + m2 

For such choice operator becomes 

^ ^ d ^ ( r"-i d \ _ 1 ^2 

8mi?V'-i ar ° V(l + r2)"-2 c»ry 2(mi + m2)i?2 " 

- \ {CDi + AD2 + 2BD3) + V{r). 



(11) 



It is obvious that 

£2 (R^ X G, K, rf) = £2 (R^, ^ £2 _ ^-^2) 

Operators Z?q, Di,D2, D3 act on the second factor in (|12|l . This action will be studied in 
the following section. 

Note that J? = for mi = m2. Let tpD £ {G,K,ii) be a common eigcnfunctions 
for operators Dq, Di, D2 if mi = m2 and also for D3 if mi 7^ m2. Then the following 
stationary Schrodinger equation 

Hifir)^PD)^Efir)ibD (13) 

is equivalent to a spectral problem for an ordinary differential equation for a function 
/(r) and an energy level E (in other words to a one-dimensional stationary Shrodinger 
equation). 2 

Proposition 2. Let be a common eigenfunction for operators Dq, Di, D2, with 
eigenvalues Sq, 61, 62 and 63 respectively. Then 

1. 5i — S2 and S3 = 0; 

2. DqiJjd is an eigenfunction for operators Dq, Di, D2, D3 with the same eigenvalues 
^0^^1,62 and S3 respectively; 

3. if DqiPjj 7^ ^pD, then DoipD^VSo'ipD oltc eigcnfunctions for operators Dq, Di, D2, D3; 
4- if DQ-tjjD ~ ipjj, then either -DoV'-D = or Si = S2 = {n — ^){n — 3)/4. 

Proof Relations [L'c-Ds] = Di - D2 and [i:'i,i:'2] = -2{Do,D3} imply 

[Do, D3]^Ijd = S3Do4'D - D3Do^D - (^1 - D2)^Ijd = {Si ~ S2)4'D, 

S3Do^Ijd + D3Do^Ijd = {Do,D3}^D = ~^[Di,D2]i^D = 0. (14) 
The last two equations lead to 

2S3Do^PD = iSl-S2)^JD- (15) 

If S3 0, then the last equation implies Dq^u ~ 4'd and the relation [Do,Di] — — 2D3 
gives S3'ipD = -DaV'i' = ^^[Do, Di]tpD — 0. Thus i^a = and equation (fT^ imphes Si — S2 
that proves the first claim of the proposition. 

Now from equation 1)14(1 one gets I?3Z?o^d — and the first two relations (jHl imply 
DiDqiPd = D2D01PD — SiDgipD- The relation DqDqi/jz) = SqDqiPu is evident, which 
completes the proof of the second claim. 

^Note that such eigcnfunctions are very special elements of the space {G,K,fi) and they do not 
span it. 
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The relation DqiI^d — SqiI^d is equivalent to (_Do + v^id) (-Do — \/5^\A) ?/'_d = 0. Now 
if DotpD 7^ V^V'-Di then t/;^ (Dq ^ \/7o id) "01) is an eigenfunction for the operator Dq. 
The function ■0^ is also an eigenfunction for operators Di , D2 , -D3 due to the second claim. 
The consideration for the function := (^Dq + \/3o id) V'n is completely similar. Thus 
the third claim is proved. 

Assume now DqiPd — Sqi/jd- Then the last relation from (jHJ gives 

It means either S'q ~ or Si ^ 62 = {n — l){n ~ 3) / 4 that proves the last claim. □ 



5 Action of operators Dq? ^1)^2,^3 in the space 

Here we use notations of sectionElfor G — SO(n + 1) and K — SO(rt — 1). Below we mean 
by the complexification of the Lie algebra g the following set 

so(n+ 1,C) = (A e 1,C)| A + A^ = .B) . (16) 

Operators Di are polynomial w.r.t. infinitesimal generators of right G-shifts. Therefore 
they conserve the spaces 7^ and generally intermix its direct summands >C^j., k — 1, . . . ,de 

with constant £ and different k. On the other hand they act in spaces Ri^i and their action 
is the same for constant £ and different i = 1, . . . , dg. 

From now we shall treat complex spaces Ri i as a simple left modules over g"'. Their 
subspaces Ri^i consist of elements annulled by the subalgebra j'' = so(7i — 1, C) C g^. 

The classification of such modules based on the notion of a dominant weight is well- 
known QI| (see also appendices A and B for a brief description). In order to apply 
this theory one should use a form of so(n + 1, C), described in appendix IXI and different 
from H16(l . Besides, since QSfc :— so{2k + 1,C) and Dk ■— so(2fc,C) are different series of 
simple complex Lie algebras, we shall consider cases of odd and even n separately. 



5.1 The case n = 2k 












J2k+1 — 





1 















In this section we shall use notations from appendix lA.il In particular, by 5Bfc we mean 
the set (jA.ip . First of all we shall construct the isomorphism g^ = 03^. in explicit form. 
Let 

i^^Ek ^^Sk\ 

e GL(2fc + l,C), 
where i is the complex unit. It is easily verified that 

J2k+lS2k+lJ2k+l = E2k+1- 

Therefore the equation A'^S2k+i + '5'2fe+iA = for A e gl(2fc + 1, C) is equivalent to the 
equation + B = 0, where B := {J2k+i) ^J^k+i- Thus the map 

B ^ 4k+iB {Jlk+iV (17) 
is the isomorphism between g"' and 
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Let 



C = 



f 





a 


A. 


a 






—a 







b 


B+ 






-B^ 










—a 


-b 




C 




V 




-Bl 






J 



where 



B+ = {bl, bk-i) , a^, 5„ a, 6 G R, C e 5o(2fc - 1). 



A- = (a_(fc_i) 



Move the second row and the second column of the matrix C to the last positions. This 
gives the matrix 



C 



/ 


A 


a A+ 


a \ 

-B^ 




-A^ 






—a 


C' 


-b 










V 


-a B_ 


b B+ 


/ 



e so(2/c + 1), C" G so(2fc- 1). 



The transformation 117|) now gives for C := J2k+i^ {■^Ik+i) expression 
/ -2ia Z_ - iZ+Sk-i V^z Z_Sk-i + iZ+ 



C ^ 



— T — T 

-Z_ - iSk-iZ. 



-V2z 



c 



V 



-Su-iZ_ + iZ^ 

Z^-iZ+Sk-i \/2z Z_S'fc_i+iZ+ 



\ 

-^/2z 
2ia / 



where Z_ := A_ + A+ + 2; ;= a + i6, C" G 25fc_i. Let us identify Lie 

algebras g"' and 53^ through the map C ^ C. Due to the definition of in section^ 
one gets the following formulas 



*12 = 


i^fefc, ^'i,fc+2 




T2 


— Fok) , '^2,k+2 






+ Fok) , 






1 

- (Fkj + J^fc,- 


-j 


+ F^kj + 


F-k-j) , j = i 


- k 


-2,3 s; 


i < fc+ 1, 






i 

- [Fkj - Fk,- 


-j 


+ F-k,- 


F-k.-j),j=i 


- k 


-2, fc + 


3 < i ^ 2A: H 


hi, 


^-2^ = 


\{F-kj+F^ 


-k, 


-j — Fkj 


- Fk-j) , j = i 


- k 


-2,3 s; 


i ^ k + 1, 




= 


- {Fkj - Fk,- 


-j 


+ F-k^-, 


- F-kj) , j = i 


- k 


-2, fc + 


3 ^ i < 2A: H 


hi, 



which imply 



fe-i 



-Di — - (i^fco - -Fofe)^ + 2 5Z + ' Fk,-] + F^k,-]} , 

1 21''"^ 
£'2 = --^ (-Ffco + -Fofc)^ + 9 X] ^F-kj - Fkj,Fk,^j - F^k,^]} , 



(18) 



fc-i 



[Fla - Fok) +^^iFkjFk-j - F^koF-k-j) , Dq = -iF, 



fcfc- 
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Since the case k = 1 does not fit the general scheme due to the triviality of the group 
K we assume from now k ^ 2. The case k = 1 will be considered below. 

Let the space TZi^i equals V^SkW for a highest weight HA.4|I . where rm £ Z+, and 
V'BkW be a subspace of V<8fc(A) annulled by the subalgebra = ^k-i- An element 
V E V'<8^(A), ^ is a highest vector of the trivial one-dimensional S^-i-module. Then 
propositions lA. ll and [^31 implv the existence of such numbers m'j G Z+, j = 1, . . . , k that 

TOfe ^ to'j, ^ TUk^i ^ . . . ^ m'2 ^ mi ^ m'l ^ —Til, 
to'j. ^ ^ m'j._i ^ ^ . . . > 77I2 ^ ^ |77i'J. 

Thus = 0, j = 1, . . . , fc — 1 and therefore mj = 0, j = 1, . . . , /c — 2. 
From now till the end of the present subsection suppose 

A = mk-iEk-i + TOfeSfc, TOfe > ruk-i > 0, mk,mk-i G Z+. 

In this case proposition I A. II implies that an every module V£i^,{m'i,ek) C VfBj.(A) contains 
the unique one-dimensional module V<s^_-^{0). This fact leads to 

dimV"<8fc(A) = mfe - r7ifc_i + 1. (19) 
Thus from proposition ^ one gets the following expansion j42| : 

£2(SO(2fc + l),SO(2A:-l),M) = {mk~mk-i + l)V<sAmkek+mk-iek-i), 

ruk ^ ruk-i 
mk,mk-i e Z+ 

where the left hand side is considered as a restriction of the left regular representation for 
the group S0(2fc -I- 1). On the other hand the space 

£2 (SO(2/c + l),SO(2fc- 1),^) 

as a Diffgo(2fe+i) (S0(2fc 4-1)/ S0(2A; — l))-modulc is expanded as 

£2(SO(2/fc-|-l),SO(2fc- 1),^) 

= (dim Vbj^ (TOfeEfe + mfe_iefc_i)) T/ojfc (m-feefe -I- mfc_iefe_i) , (20) 

rUk ^ TOfc-i 
mk,mk-i e Z+ 

where dim V®^ (rnkSk + mk-iSk-i) is given by (|A.8p . 
Let 

fe-i ^ fc-i ^ 

fc-i 

C := C'l£2(so(2fc+i),so(2fc-i),p) = ^fefe + {-f^feo, -Fofe} + ^ ({-f^fcj, i^jfc} + {-Ffc -j, F_jfe}) 

i=i 

be operators from Diffgo(2fc+i) (S0(2fc + 1)/ S0(2fc — 1)), where C is the universal Casimir 
operator ljA.6|) . Due to (|A.3p and l|A.9|l the operator £>+ "raises" weight subspaces of 
V(8j^(A) and the operator D~ "lowers" them. 

Since [FkjTFk~j\ — [F^kjiF^k-j] — one gets the following relations 

D^^D++D- + ]^ (fI, -C),D2^-D+ -D- + - c) , Z?3 = i - D') , 
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D^^l (Di - D2) - D- {Di - D2) + C = -Dl -D,- D2. (21) 



Commutator relations (jHJ now give 



[Fkk,D+] = 2D+, [Fkk,D-] = -2D-, (22) 
[D+,D-] = -^F^^ + ^CFkk + \{2k-l){2k-3)Fkk. (23) 



Formulas ljA.5|) and (|A.7p implies 
C 



iV /, /, 1\^ /, 3^^^ 



i>.,(.) yy' + ^^-2) +^fc + --i-2j -['^-2) -['^-2) jid- (24) 

It follows from the paper PH) that 

Vs, (A) = ® l^-(.-2)e, ® ■ • ■ ® 1^(.-2)s. ® Ke. , (25) 

where v — — rrik-i and all summands are one-dimensional weight spaces w.r.t. the 
Cartan subalgebra Ifk-^ Formulas (|A.3p and ljA.91) imply 



The action of operators Fkk, , in the space V<Sk (A) was calculated in w.r.t. 
some base. In particular, in Vts^W there are no nontrivial invariant subspaces w.r.t. this 
action. We shall obtain simpler formulas for the and D^-action w.r.t. a base in Vqjj. (A) 
with a normalization different from those in |3f)j . 

Lemma 1. Let := {—i>, —i> + 2, . . . , v — 2, v). There is a base {xj)j^j^ inV<s^.{X) such 
that 

FkkXj = JXj, D+Xj = ^(i ~rnk~ mfe_i - 2k + 3){j - v)xj+2, (26) 
D-x, = \{]+ + ruk-i +2k- 3)(j + i^)xj-2, (27) 
where Xj — if j ^ Li,. 

Proof. Since the action of an algebra, generated by operators Fkk, , , is irreducible 
in V<8^(A), one can define by induction nonzero elements Xj € V^e^, j € Li, such that 
formulas 1)26(1 are valid. Prove by induction formula (|27|) . For j = —v it is evident. 
Suppose that H27() is valid for j = —v, — + 2, . . . , i, where i <v. Then using (|24|l one gets 

i(i - mfe - mk-i -2k + 3){i - y)D-x^+2 = D-D+x^ = {[D-,D+] + D+ D') 

\Flk - \cFkk - \{2k - \){2k - 3)Fkk^ x^ 
+ -{i + mk+ nik-i + 2k- 3){i + v)D^Xi-2 

(""^ + ""'-1 + " + " i)mk-i + ^(2^ - l)(2fc - 3)^ ^ x^ 
+ ^(j + mfc + r7ifc_i + 2fc - 3)(z + v)[i - ruk - rrife-i - 2fc + l)(i - 2 - j/)xi 



^In appendix |D] we shall give a proof of expansion 1251 independent from the theory of Yangians used 
1^ . 
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= - ruk- TUk-i -2k + 3){i - + nik + nik-i + 2k - l)(i + 2 + u)xi, 
due to the identity 

(i — rrife — TUk-i — 2k + — h'){i + rrik + rrik-i + 2fc — + 2 + 1/) 
— {i + rrik + rrik-i + 2k — 3)(2 + — nik — rrik-i — 2k + — 2 — v) 

Since (z — — ruk-i — 2k + — i^) 7^ we obtain 

= ^(i + "i/c + mk-i +2k-l)(i + 2 + v)xt 

that completes the induction. □ 

Lemma ^ expansion (|20|l and relations (|21|l effectively describe the action of oper- 
ators Dq, Di, D2, D3 in the space £^ (S0(2fc + 1), S0(2fc — 1), ^i)- Consider the problem 
of finding all common eigenvectors ipD of operators Dq, Di, D2 and optionally D3. It is 
equivalent to the problem of finding all common eigenvectors of operators Z?g, + 
and optionally £>+ — in the space V<s^{X)- 

Eigenvectors for the operator Dq are 

c+Xj + c~x-3, c± e C, j e L^, j ^ 

with eigenvalues —j^. Since 
(D+ + D~) (c+Xj + c-x-j) = ^(j -ruk- mfc_i - 2fc + 3)(j - j/) (c+Xj+2 + c-X-1-2) 

+ ^{j + mk+ mk-i + 2k- + v) {c+Xj~2 + C-X^]+2) , 

the requirement 

+ D-) {c+Xj + c-x-j) - c+xj + c-X-j 

implies {j — mk — rrik-i — 2fc + 3)(j — v) = that leads to two cases: j ~ — nik-i and 
j = rrik + ruk-i + 2fc — 3. 
In the first case one gets 

= (TOfe - mfc_i) (^TTlk + ^ - (c+Xm;,-mfc_i-2 + C-X- 

that implies one of three possibilities 

1. ruk - ruk-i = 0; 

2. rrifc - mk-i - 2 = -rUk + mk-i] 

3. mfc — mk-i — 2 = 0, c+ + c_ = 0. 
Thus we obtain the following eigenvectors: 

1. {D+ + D-) xo = for nik ~ ruk-i = 0; 

2. (1)+ + D-) (xi + X-i) ^ {mk + - i) (Xi + X-i) for ruk G N, rrik-i = nik - I; 
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3. {D+ + D ) (xi - x-i) = - (m-fe + - I) (Xi - X-i) for "ifc e N, mfc_i = Wfe - 1; 

4. + D-) (x2 - X-2) = 0, ruk-i - mfe - 2, = 2, 3, . . .. 

In the second case one gets ruk + rrik-i + 2k — S = j ^ irik — rrik-i that implies 
^ rrik^i ^ I — A; and thus k = 1 that contradicts to the assumption k ^ 2. 

Using relations (|21|l this consideration can be summarized in the following proposition. 

Proposition 3. For n — 2k, k ^ 2 there are four series of common eigenvectors in 
VsJ^fcEfc + mk-iSk-i), mfe,mfe_i G Z+ for the operators Dl,Di,D2: 

1- D^xo = D3X0 = 0, Dixo = D2X0 = -rrikimk + 2k ~ 2)xo, mk = rrifc-i; 

2- Dlixi + X-i) - -ixi + X-i), D2U1 + X-i) = ~mk{mk + 2k ~ 2)(xi + X-i), 
Diixi + X-i) = {-ml - 2{k - 2)mk + 2fc - 3) (xi + X-i), 

^3(xi + X-i) = i {mk + - I) (Xi - X-i), mk-i = ruk - 1, ruk G N 

3. Dlixi - x-i) - -ixi - X-i), Diixi - x-i) = -mkim,k + 2k - 2)(xi - X-i), 
^2(xi - = {-ml - 2{k - 2)mk + 2fc - 3) (xi - X-i), 
Dsixi - X-i) = -i {mk + - i) (Xi + X-i), mk-i = nik - 1, G N; 

4- Dl(x2 ~ X-2) = -4(x2 - X-2), £'3(X2 - X-2) = -4i (mfc + /c - §) xo, 

^i(X2 - X-2) = D2{X2 - X-2) = {-ml - 2{k - 2)mk + 2/c - 3) (x2 - X-2), 
nik-i = mfe - 2, mfc = 2, 3, 4, . . . 

Only the first vector is also an eigenvector for the operator D^. 

Multiplicities of corresponding eigenvalues in (SO(n + 1), SO(ri — 1), //) are equal to 
dimV<8fc {mkSk + ruk-iSk-i) and can be calculated in explicit form using fA.tS)) . 

Consider the case fc = l,n = 2. Now the group K is trivial and therefore Vf8i(A) = 
V<8i(A). The algebra *Bi — so(3,C) ~ s[(2,C) = 2li is spanned by elements Fii,Foi,Fio 
with commutator relations 

[Fii,Foi] = —Fqi, [Fii,Fio] = Fio, [-Fic-Foi] = -^ii- 

Its representation theory is well known: all its finite dimensional irreducible modules are 
of the form 

where m S Z+ U (Z+ + i), aU Vj^^ are one-dimensional weight subspaces w.r.t. [}i = 
span(Fii) and 

^^10 : Vje^ > j = -m, . . . , m - 1, Foi : Vje, , j ^ -m + 1, . . . ,m 

are bijections. 

We shall consider only m G Z_|_ since 

£2(SO(3),Ai) = (2m + l)%,(m£i). 

Thus there are additional weight subspaces in the module Vfs^{mei) w.r.t. expansion (|25(l 
and the action of the algebra, generated by the operators -D+ = ^F^q, = ^Fq^, is not 
irreducible in Vts-^{mei). 

One can choose a base (Xi)JL„m Vts-^ijnei) such that 



Xj e Vjei, F11X3 =jXj, FwX] = ^ V(™ - i){m + i + l)Xi+i, 
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FoiXi = --/^V{m + j){m-j + l)Xj-i, 

where as above Xj = for |j| > m. 

Eigenvectors for the operator Dq = — -fn are 

c+Xj + c-X-j, c± e C, j = 0, 1, . . . ,m 

with eigenvalues — j^. Since 

{D++D-) (c+xj+c^x-j) 

= ^V(m - j)(m + j + j - l)(m + j + 2) (c+Xj+2 + c_X-j-2) 

+ ^ V(™ + J)("^' - i + + - - J + 2) (c+Xj-2 + C-X-j+2) , 

the requirement 

(£>+ + D") (c+Xj + c^X-j) ~ c+Xj + c-x-j 
impUcs (m — j){rn + j + l){m — J — l)(rn + j + 2) = that gives two cases: j = m and 

J = TO — 1 . 

In the first case one gets 

(£>+ + D^) (c+Xm + C-X-m) = ^ V n^(2TO - 1) (c+Xm-2 + C-Xm+2) 

that impUes one of three possibilities 

1. TO = J = 0; 

2. TO — 2 = — m; 

3. TO - 2 = 0, c+ +c_ = 0. 

This gives the following eigenvectors: 

1. {D+ + D-)xn = 0, TO = 0; 

2. + D-) (xi + X-i) = 1 (Xi + X-i) , m = 1; 

3. + D-) (xi - X-i) = (Xi - X-i) , m = 1; 

4. (I?++i?-)(x2-X-2) =0, TO = 2. 

It is easily seen that these eigenvectors corresponds to eigenvectors from proposition O 
for nik = TO,TOfc_i = 0. 

In the second case it holds 

{D+ + D^) (c+Xm-l + C-X-m+l) = ^ V3(2to - 1)(to - 1) (c+Xm-3 + C-Xm+s) 

that implies one of three possibilities 

1. TO = 1, j = 0; 

2. TO — 3 = — m + 1; 

3. TO — 3 = 0, c+ + c_ = 0. 

Thus one gets the following eigenvectors: 
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1. iD+ + D-)xo = 0, m=l; 

2. {D+ + D-) ixi + x-i) = 1 (Xi + X-i) , m = 2; 

3. {D+ + D-) (xi - X-i) = -| (Xi - X-i) , m = 2; 

4. (I?+ + i?-)(x2-X-2) = 0, m = 3. 



Since it holds C 



= m{m + 1) id and relations 12111 are valid also in the case 

Kg J (me 1) 

k = 1 one gets the following proposition. 



Proposition 4. There are eight common eigenvectors in V<Si{rnei) for the operators 
DlDi,D2: 



1. 


Dlxo-- 


= Dixo 


= -D2XO = £'3X0 = 0, TO = 


0; 


2. 


Dlxo-- 


= D3X0 


= 0, DiXo = £'2X0 = -Xo 


TO 


3. 


Dlixi 


+ X-i) 


= f2(xi +X-i) = -(xi + 


X- 




Dsixi 


+ X-i) 


= ^(Xi - X-i), m^l; 




I 


Dlixi 


-X-i) 


= £'i(Xi - X-i) = -(xi - 


X- 




Dsixi 


-X-i) 


= -^(Xi +X-i), TO = 1; 




5. 




- X-2) 


= -4(X2 - X-2), Di{x2 - 


X-: 




Di{x2 


- X-2) 


= -\/6ixo, TO = 2; 




6. 


Dlixi 


+ X-i) 


= ^i(xi +X-i) = -(xi + 


X- 




Dsixi 


+ X-i) 


= ii(xi - X-i), TO = 2; 




7. 


Dlixi 


-X-i) 


= £'2(xi - X-i) = -(xi - 


X- 




Dsixi 


-X-i) 


= -fHxi +X-i), m = 2; 




8. 


DI{X2 


- x-2) 


= £'l(X2 - X-2) = -D2(X2 


- X 






- x-2) 


~ — \/30ixo, TO = 3. 





-(X2 - X-2), 



-4(xi-X-i), 



Only the first and the second vectors are also eigenvectors for the operator D3. 
Multiplicities of corresponding eigenvalues in (S0(3),/i) are 2m + 1. 

5.2 The case n = 2A; - 1 

Here we use notations from appendix lA. 21 The algebra Sfc is considered there as a subal- 
gebra of *Bfc. Therefore one can easily obtain analogs of formulas H18|l simply by deleting 
the terms Fko and Fok'- 

^ k-l 

i=i 



■^2 = 2 X! ^^-^J ~ Fkj,Fk-j - F^k-]}, D3 = i^ (Ffcjffc _j - F^kjF^k-j) ■ 
i=i i=i 

Let the space TZi^i equals Vs)^(A) for a highest weight (|A.10|) . where mi e Z+, i ^ 
2,TOi e Z, and Vi)^(A) be a subspace of Vs^W annulled by the subalgebra fi"- = Dk-i- 
Reasoning as above in the case n — 2k, one gets that V^2)j.(A) is nontrivial iff 

A = TOfc£fe + TOfc_i£fc_i, TOfc > |TOfc_i|, TOfc G Z+ , TOfc_i G Z^, (28) 
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where Z'j. = Z+ for fc ^ 3 and Z2 = Z. In this case one has 

dimys,(A)=mfe-|mfe_i| + l. (29) 

Below in the present subsection we suppose that condition (|28|l is vahd. This leads to 
the expansion: 

£2 (SO(2fc),SO(2fc-2),^) 

= ^ {mk-\mk-i\ + l)Vz)^{mkek + mk-iek-i) 

rrik ^ |mfc_i| 
ruk e Z+, m/c-i e ZJ. 

of the left S0(2A:)-space (S0(2A:), S0(2A: - 2), /i) and to the expansion: 

/:2(SO(2fc),SO(2fc-2),^) (30) 

= ^ {dimV's^ {nikek + mk^iek-i))Vs^ {mkSk + mk-iSk-i) , 

nik ^ |mfe_i| 
mk , ruk-i e 1,'k 

of the same space as a Diffgo(2fc)(SO(2fc)/ S0(2fc — 2))-module, where the dimension 
dimVj)^ (mfcEfc + mk-iSk-i) is given by (|A.8(I. 
Now let 

fc-i fe-i 

C := C'|£2(go(2fc),so(2fc-2),p) = -^fefc + ( { -^^i ' i + {-^fe -j'-^-ife}) 

i=i 

be operators from Diffso(2fc) (S0(2fc)/ S0(2fc — 2)), where C is the universal Casimir op- 
erator ljA.ll|) . 

Formulas (EU and ^ are valid without any modification and formula (|23() becomes 
[D+,D-] = -iit, + icFfefc + (fc - l)(/c - 2)Ffefc. 



Now 



C 



{{nik + k-lf + {mk-i +k-2f -{k- if - (fc - 2)^) id . 



From it follows^ that 

Vs, JA) - ® 1^-(.-2)s. ® ■ ■ • ® V'(.-2)e, ® Ke. , (31) 

where v = rrik — \mk~i\, all summands are one-dimensional weight spaces w.r.t. the Cartan 
subalgebra l)fc C Dk and the algebra, generated by the operators D~^,D~, acts in Vs)j^(A) 
in an irreducible way. 

Again we shall simplify formulas for this action w.r.t. 35^ using another base. The 
next lemma can be proved completely similar to the proof of lemma ^ 

^See also appendix [H] for a proof independent from i35l . 
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Lemma 2. Let v := rrik — \mk-i\, Ly := {—v, —v + 2, 2, i/). There is a base 
iXj)j(zL^ V^^{X) such that 

FkkXj = jXj, D+X] = ^{i -fnk- \mk-i\ - 2k + 4)(j - i^)Xj+2, 

D^Xj = ^(i + ruk + \mk-i\ + 2fc - 4)(j + i^)Xj-2, 

where Xj — if j ^ Ljj. 

Arguing as in the *Bfc-case one gets the fohowing proposition. 

Proposition 5. For n — 2k — 1, k ^ 2 there are four series of common eigenvectors in 
Vs^imkEk + \mk-i\sk-i), nik G Z+, mk-i G Z'j, for the operators Dl,Di,D2: 

1- DIxq = -D3X0 = 0, DiXq = D2X0 = -mk{mk + 2k - 3)xo, mk = \mk-i\; 

2. Dlixi + x-i) = -ixi + X-i), D2{xi + X-i) - -mk(mk + 2k - 3)(xi + X-i), 
Diixi + X~i) = {-rnl + (5 - 2k)mk + 2k - 4} (xi + X-i), 
Dsixi + X-i) = i {fnk + fc - 2) (xi - X-i), = "ifc - l,mfe e N 

3- -Do(Xi - X-i) = -(Xi - X-i), -Di(xi - X-i) = ~mk{mk + 2k - 3)(xi - X-i), 
-D2(xi - X-i) = {-ml + (5 - 2fc)TOfe + 2fc - 4) (xi - X-i), 

-D3(xi - X-i) = -i(™fe + fc - 2) (xi + X-i); \fnk-i\ = mk ~ l,mk £ N; 

4- Dl{x2 - X-2) = -4(X2 - X-2), £'3(X2 - X-2) = -41 {mk + k- 2) xo, 
i^i(X2 - X-2) - D2{X2 - X-2) = {-ml + (5 - 2fc)mfe + 2fc - 4) (x2 - X-2), 
|mfc_i| = mfc - 2,mfc = 2,3,4, .. . 

OnZ?/ the first vector is also an eigenvector for the operator D^. 

Multiplicities of corresponding eigenvalues in (SO(n + 1), SO(ri — l),/i) are equal to 
dim Vx)^ (rakEk + mk-\Ek-\) and can be calculated in explicit form using iA.S\) . 

Remark 1. For k ~ 2 a value rUk-i = mi can has an arbitrary sign and one gets eight 
common eigenvectors found in \S4Jl - 

Remark 2. Results of propositions\^^and\^ correspond to proposition\^ and are even 
more restrictive. Indeed, if ipn G £^ (SO(n + 1), SO(n — 1), /i) is an eigenf unction for 
operators Dq, Di, D2 and D^, then DqiPd — D^ipo — 0, -DiV"!) — D2'(pD- 

6 Scalar spectral equations and some energy levels for 
the two-body problem in S" 

Here we shall consider the spectral problem (|13|l . where the operator H is defined in 
and V'-D is one of common eigenf unctions for operators Dq, Di, D2 and optionally D^. In 
this section m denotes reduced mass HlOfl and integers ruk correspond to highest weights 
in 5o(n)-modules. 

Let DqiPd = SoipD, DiipD — Siipu, i = 1,2. In accordance with remark[21there are two 
main cases: 

1. D^-tpD = 0, Sq = 0, Si = 62, particle masses are arbitrary; 

2. D^ipD 7^ ipu, particle masses are equal. 
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In the first case 



{CDi + AD2 + 2BD3) =h{C + A)^D 




In the second case 



^ = 1 — i( = U, C = - — 
In all cases one gets the following spectral equation for the function f{r) 




(32) 



a, 6, c 0, < r < 00. 

where coefficients a,b,c are described below. 

For eigenfunctions ipn classified in proposition {n = 2k, fc = 2, 3, . . .) one has 

1. a = c = mk{mk + 2fc — 2)/8, b = 2a, G Z-|_, masses are arbitrary; 

2. a = TOfe(TOfe + 2/c - 2)/8, b = {ml + {2k - 3)mk - fc + 2)/4, c = (m| + 2(fc - 2)mfc - 
2fc + 3)/8, TOfe G N, masses are equal; 

3. a = {ml + 2(fc - 2)mfc - 2fc + 3)/8, b^{ml + {2k - 3)mk - fc + 2)/4, c = mk{mk + 
2k — 2)/8, mfe e N, masses are equal; 

4. a = c = (m| + 2(fc-2)mfc-2fc+3)/8, b = {ml + 2{k-2)mk-2k+5)/4, = 2,3, . . ., 
masses are equal. 

Proposition^ {n — 2) gives the following values for a, b, c: 

1. a = c = b = 0, masses are arbitrary; 

2. a = c= 1/8, 6 = 1/4, masses are arbitrary; 

3. a = 1/8, 5 = 1/4, c = 0, masses are equal; 

4. a = 0, 5 = 1/4, c = 1/8, masses are equal; 

5. a = c = 1/8, b = 3/4, masses are equal; 

6. a = 1/2, b = 3/4, c = 1/8, masses are equal; 

7. a — 1/8, b = 3/4, c = 1/2, masses are equal; 

8. a = c = 1/2, b = 3/2, masses are equal. 

Finally, proposition |S1 corresponds to the following cases (rt = 2fc — 1, fc = 2, 3, . . .) 

1. a = c = mk{mk + 2fc — 3)/8, 6 = 2a, G Z_|_, masses are arbitrary; 

5 

2. a = mfc (rrifc + 2fc- 3)/8, b = (m^ + (2fc -4)mfe - fc + -)/4, c = (m^. + (2fc - 5)TOfc - 
2fc + 4)/8, TOfc G N, masses are equal; 

5 

3. a = {ml + {2k - 5)mk ~2k + 4)/8, b = {ml + {2k - 4)rnfe - fc + -)/4, c = mk{mk + 
2fc — 3)/8, mk G N, masses are equal; 

4. a = c= (m2 + (2fc-5)mfc-2fc+4)/8, 6= (m2. + (2fc-5)TOfe-2fc + 6)/4, = 2,3,..., 
masses are equal. 

We shall consider equation for the Coulomb and oscillator potentials. 
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6.1 Coulomb potential 

For the Coulomb potential 



14 = -^cot-^ = (r- i ) , 7 > (33) 



R R 2R\ r 

theorems [T1 and IB. II imply the self-adjointness of the two-body Hamiltonian Hy^ with its 
domain defined by (|B.Hl . where l/i = for < r < 1 and Vi = for 1 ^ r < oo. 

Equation l|32(l for V — Vc is the Fuchsian differential equation (see appendix Q with 
four singular points r ~ 0, ±i, oo and corresponding characteristic exponents: 

P± = ^ (2 -n± V("-2)2 + 32a) , p^^'' = i (2 ~ n ± ^ (n - 2)^ + 32c) , 



^ - ( n - 1 ± y (n - 1) + 8 (TO£'i?2 _ hnR-/ + a - b + c) ] , (34) 



= i - 1 ± ^{n-lf +8 {raER'^ + imR-y + a - b + c) 

Here and below we suppose that a square root for a positive number is positive; for other 
numbers it is an arbitrary root. 

The requirement f{r)^D G T>OTa{HvJ) restricts asymptotics of /(r) near singular 
points r = and r = 00. Let /(r) ^ r'''"' as r ^ +0 and /(r) ^ r^''*°°' as r ^ +00. We 
shaU show that /(r)V'n G Dom(i7vJ iff pf") = and = 

The inclusion 



evidently implies p^°^ > —n/2, p^°°'> > —n/2. On the other hand one can easily see that the 
inequality a ^ 1/8 leads to p^^' ^ —n/2 and the inequality c ^ 1/8 leads to p^'' ^ —n/2. 

From the consideration above it follows that if a < 1/8 (c < 1/8) then a = (c = 0). 
For a — the inequality p^^^ = 2 — n > —n/2 implies n < 4. 

(0) 

For a = 0, n = 3 the asymptotic /(r) ~ r''- — 1/r means that A(/i/'/)) ~ S{0) as 
r ^ that contradicts to 

A(/V.i3)e£L(S"xS",xxx), (35) 

see theorem lB.il 

For the case a = 0, n = 2 it holds p^_^'^ = = and the theory of Fuchsian differential 
equations jSlJ, [Ml implies that canonical asymptotics of a solution for (|32|l near r = 
are 1 and logr. The latter asymptotic again leads to A(/7/)/)) ^ 5(0) as r ^ that again 
contradicts to IpT^ . 

(0) 

Thus in all cases it should be /(r) ~ r''+ as r ^ 0. Reasoning in the similar way one 

_ (00) 

gets also in all cases the asymptotic f{r)^r as r ^ +00. 

Consider the problem of reducing equation H32|l with potential H33(l to the hypergeo- 
metric equation via reducing (|32|l to the Heun equation by transformations (|C.2|I . lC.3p 
and then using thcorcm lC.il 

Singular points of equation 13211 form a harmonic quadruple (see appendix Q. There- 
fore, one can use only the first case of theorem lC.il Move singular points (0,±i,cxD) of 
equation H32|l to the quadruple (0, l,2,(xi) by a fractional linear transformation t — rlr) 
of independent variable. 

Since the order of singular points on a circle or on a line is conserved by such trans- 
formation only two possibilities can occur. The first one corresponds to the map of the 
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unordered pair (±i) into the unordered pair (0,2). The second one corresponds to the 
map of the unordered pair (0,oo) into the unordered pair (0,2). 

Then one can reduce the transformed equation to the Heun one by a substitution of 
the form ljC.3(l . One of requirements of the first case of theorem IC.ll is the equahty of 
characteristic exponents at points and 2. In terms of characteristic exponents HH4|I it 

means that either jp^^ - = \p'+''^ - pL"'^! or \p^°^ - = 1/9^°°^ " P'-\- The first 
possibihty can not occur for a nontrivial 7. Therefore, not loosing generahty, one can 
consider the map 

t = T(r):=^, r: (-i, 0, i, 00) ^ (cx), 0, 1, 2). 
r + 1 

This map transforms equation (|32|l with potential H33|) into the equation 

ht{t)+A{t)!tit)-B(t)!(t) = ^, |i-l| = l, lmt<0, (36) 

where 

nt^ -2nt + 2n-2 



A{t) 

B{t) = 2 



t{t-l){t-2) ' 

(ERH^t - 2)2 + i?7it(t - 2)(t2 -2t + 2)) + a{t - 2)* - bt\t - 2f + ct^ 



i2(t- l)2(i_2)2 
The substitution 

transforms (|36(l into Hcun equation (|(].10(l with the parameter 7' instead of 7, where 

a = pf + p^^ + pi-^ + p^-^ p = pf + pf + + pi-^ d = 2, 
7' = l-pL°'+pf ,<5 = l-p«+p«,£ = l-p(!-)+pf). 

(0) (i) (00) 

Here t''+ {t — \Y+ [t — 2Y+ means the function holomorphic on C\(— 00, 2] and real for 
real t > 2. Restrictions on asymptotics of the function / near the points r = 0,oo are 
equivalent to the boundedness of the function w{t) near the points t = 0,2. 
Obviously, the accessory parameter q can be found as 

g = -21imi -S(t)+ ^ + /±- + ^L4(i) ' ' 



t t-l t-2 j ^ ' i2 t{t-l) 

r, (0) (00) \ 

+ j = ^Pfp'^ + 2pfpr -in- 3)pf + (n - l)(2p« + p^^^) 

- 4mRYi + 16a. (37) 

Theorem IC.ll implies that this Heun equation can be transformed into the hypergeo- 

metric equation by a rational change of independent variable t z : z = P(r), where P 
is a rational function, iff 

y = e, (38) 

aP-q = 0. (39) 

Equation 1)38(1 is equivalent to 

a = c. (40) 
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Using the equalities 

a = pf^ + p^i^ + + i - 1 + ^{n-lf + 8 {mEK^ + imR-i + a - b + c)^ , 



(3 = + + + 1 _ 1 _ {n -1)2 + 8 {mER^ + imi?7 + a - 6 + c) j , 
one can rewrite equation as 

(^pV"^ + p^' + + i(n - 1)^ - i ((n -1)2 + 8 (m£;i?2 + imi?7 + a - 6 + c)) 

- 4pV"V^^ - 2pf + - 3)pV"' - (« - 1)(2P^^ + + 4mi?7i - 16a ^^^^ 

(p«)% fpi-')' + 2p« (pV-) - pf ) + (2n - 4) pf 



- (ra - 1) p^' + 2mi?7i - 2mER^ - 18a + 26 - 2c = 0. 
Excluding squares of values p^' , p^' , p^'' from \A1^ with the help of obvious equations 

[pf)\{n-2)pf-Sa = 0, 
(p^^) -{n- l)p^' - 2mi? (i?^; - 7i) - 2(a - 6 + c) = 0, 
(pV°°')' + (n-2)p(r^-8c^0 
for characteristic exponents, one gets 

(2p« - n + 2) (p^f - pf ) + 8(c - a) = 0. 

For a = c it holds p^"* = p+'' and thus equation (|39|l is a consequence of 140|l . 

From now till the end of the present subsection we suppose that a = c. This condition 
corresponds to cases 1,4 of proposition cases 1,2,5,8 of proposition 0] and cases 1,4 of 
proposition [S] 

The fist case of theorem lC.ll implies then that the function w w.r.t. a new independent 
variable 

z := 1- (t- 1)2 ^i(2-i) (42) 
satisfies the hypergeometric equation: 

z(l - z)w"{z) + (7 - (5 + ^+ l)z)w'{z) - 5/3w(z) = 0. (43) 

The correspondence between characteristic exponents of the Heun and the hypergeo- 
metric equations connected by H42|l implies 

1 1 1 n :tt-„ TT- 1 



7 = 7' = 1 + V(^i-2)2 + 32a G M, 5 = -a - + - ^(^i - 2)^ + 32a + -(s + s) G R, 



where s = y (n - 1)^ + 8 (jnER'^ + imi?7 + 2a - 6). 
Since 

.\ 2 

, r — I 

z-l 



r + i 
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the half-line [0, oo] on the r-planc is mapped into the circumference on the z-plane defined 
by the equation |z — 1| = 1, while the values r = 0, cxd correspond to the point z — 0. 

The function w{z) is bounded near the point z and 1~7 = -^(n - 2)2 + 32a < 0; 
therefore it holds (see l|C.5|l ^ 

wiz) = w^{z) := c+F(5, /3; 7; z), z e {\z — 1| = 1, Im z > 0) , c+ = const, 
w{z) = w- (z) := c_F(S, /3; 7; z), z G (|z — 1| — 1, Im z < 0) , c_ = const . 

An equivalent problem for the hypergeometric equation was considered in |14|. However 
there was made an assumption equivalent to c+ — c_ without any proof (formula (3) in 
|14j). Below we fill this gap. 

Functions w±(z) should be analytic continuations of each other through the regular 
point z = 2.^ Due to formula (jC.6p (applicable since j — a — (3 ^ R) it means that functions 

rf^)rf^ -a- 3) 

c+ ^l' — W^F(5,/3;a + /3-7 + 1;1 - z), |z - 1| = 1, Imz > 

r(7 - a)r(7 - /?) 

and _ 

^_^dLAl ^F(5,/3;S + /3-7 + l;l-z), |z - 1| = 1, Im z < 

r(7-5)r(7-/?) 

are analytic continuations of each other through the point z = 2 as well as functions 

c+^-ffi^^-t£— ^(l-z)^-"-^F(7-5,7-/3;7-5-/?+l;l-z), |z-l| = 1, Imz >0 
r(5)r(/3) 

and 

c- ^ ~ (1 z)^-"-^F(7 - a, 7 - /3; 7 - a - /3 + 1; 1 - z), |z-l| = 1, Imz < 0. 

r(5)r(/3) 

The first requirement is equivalent to the equality 

(c.-c_)lffl:fc^=0, (44) 

r(7-5)r(7-/3) 



while the second one to the equality 



'\ " "^'JJ r(5)r(/3) 

Since 7 — a — /? ^ M linear system H44II , (|45|l has a nontrivial solution , c_ iff either 

r(7)r(7 -5-/3) r(7)r(5 + ^-7) „ 

— = or — ~ 0. 

r(7-5)r(7-/3) r(5)r(/3) 

Taking into account 7 — /3 ^ M, /3 ^ M, one gets 7 — 5 = — fc + 1 or 5 = —k +1, fc e N. 

Not loosing generality suppose that Re s < 0. Then the first equality is impossible and 
the second one yields 

4im_R7 



s = 1 - 2fc - ^/{n - 2)2 + 32a + 



l-2k- ^(n- 2)2 + 32a' 



^Recall that the function -F(a', /?'; 7'; z) is holomorphic in C\[l, +00). 
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since Im — SimR'j. From the definition of s one gets therefore the following formula for 
energy levels: 

= 37^2 ( o (fc' - + 1) - T + 2a + ^ + —r- Vin - 2)2 + 32a j 



TOi?2 V2' 'A 4 

2m72 



(v/(n - 2)2 + 32a + 2fc - l) 



These energy levels are degenerated and their multiplicities coincide with multiplicities 
of eigenvalues in propositions 01 El and El 

Taking into account all transformations used while reducing equation H32|) to the hyper- 
geometric one we get the following expression for radial eigenfunctions (up to an arbitrary 
constant nonzero factor) 



+ i)Pr+p';'+pV~' ^ - J - 1)! (7), {r + i)2^" ' 



where p'^\ pl^\(3 and 7 are given by above formulas for E = Ek- 
6.2 Oscillator potential 

The oscillator potential for the sphere S" has the form 



K(r) = -RW tan^ ^ ^ J^^^' ^ ^ 



It has a positive singularity along the sphere equator and looks like an infinite potential 
well. Therefore from the physical point of view it is natural to consider wave functions 
defined on M' and vanishing as r — + 1. 

From the mathematical point of view theorem IB. II is not applicable since 



yo^£L(S"xS",xxx). 

use the Friedrichs extensior 
the domain given by theorem IB. 21 where M' C S" x S" is defined by the inequality 
r = tan(p/(2i?)) < 1. 

Equation H32|l for y = Vo is a Fuchsian one with six singular points O,±l,±i,oo and 
corresponding characteristic exponents: 



P± - o (2 - n ± V(^-2)2 + 32aj , p^^^ ^-{2-n± ^ {n - 2)2 + 32c j , 







pT 


= pi-' 


pi' 


- pi-" 



1 



pi' = P±'' = ^ (1 ± V^l + 4i?4mu;2) _ 



Similarly to the previous section the function /(r), r € (0,1) should be ~ r''+ as 
r +0. On the other hand the inclusion 

implies the convergence of the integral 

92 (V(7^, V ifijD)) dx X dx (46) 
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where (72 is defined in jSJl and V means the gradient operator. 

The convergence of H46|) is equivalent to the convergence of its " radial part" 



(1 + r2)"-4- 

Therefore if / r'''^' as r 1 - 0, then p^^^ > 1/2 and thus p^^^ = r''+\ 

Conversely, it can be easily verified that if / is a solution of %VZ\ with asymptotics 

(0) (1) 
/(r) ~ r''+ as r — > and f{r) ^ r''+ as r — > 1 — for V ~ Vo, then f{r)'iljD S 

Dom((i7yJ^). 

Fortunately, one can glue points r — ±1 together (as well as points r = ±i) by the 
change of the independent variable r — > ^ = r^, which transforms the differential 
equation under consideration into the following Fuchsian differential equation with four 
singular points: 

2C(C + i) •^^ + C(CTTFr^ r~(C^j~C-'"^9^^'' (47) 
< c < 1- 

Singular points — 1,0,1, c» of this equation form a harmonic quadruple and correspond 
respectively to characteristic exponents: 

Ji) 1(0) (1) 1 

P± ' ' P± ' 2 ± ' 

The same arguments as for the Coulomb problem leads to the conclusion that the only 
possibility to transform equation H47|l to the hypergeometric one via transformations (|C.2|I . 
(|C.3|I and then using theorem lC.il corresponds to the map of the unordered pair (0,cxd) 
into the unordered pair (O, 2) by a Mobius transformation. 
Not loosing generality, one can consider the substitution 

^ = ^(0 = ^, ^: (-1,0,1, cx.)^ (00,0,1,2). (48) 

The interval under consideration for the variable t is again (0, 1). Substitution (|48|l trans- 
forms equation (|47ll into equation H36() with 

n(i-l) , 2 / 0/ R^uj'^t{t-2)\ 2a , ct 

Define a function w{t) by 

1 (0) (1) 1 

It satisfies Heun equation (|C.10() . where 

a = \pf + P« + Ipi-' +f^,f3^ IpT + + Ip'r' + = 2, 

1^(0) Jo)\ x_i , Ji) Ji) I 1/^°°) 



7=^+2 V+ ^ P- ) ' = ^ + Z'- ' ^ = + 2 V + P- 

_1 (0) _ (1) _1 (00) 

Here i 2P+ (t — 1) (i — 2) 2P+ means the function holomorphic on the domain 
C\((— oo,G]U [l,+oo)) and real for real t G (0,1). Restrictions on asymptotics of the 
function /(r) near the points r = 0, 1 are equivalent to the boundcdness of the function 
w{t) near the points t = 0, 1. 
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Calculation, similar to (|37|l . yields the following value of accessory parameter q for 

g = -2mR^E + 2b + n + Ip^^^) + p'^' + ^pT pi°"' + jpT ■ 

Condition 1)38(1 of theorem IC.ll is again equivalent to 1)40(1 . Condition 1(39(1 of the same 
theorem can be written as 

which is again a consequence of 1(38(1 . 

Suppose that condition ((38(1 is valid. Thus we are in the situation of the first case of 
theorem IC.ll and changing the independent variable t by a new one z according to 1(42(1 . 
one gets hypergeometric equation ((I^ with 



5 = -a = - (^2 + ^(^^"2)2+^ + Vl + 4i?4TO^2 ^ 5 j ^ 
P^\P= \ (2 + \/{n - 2)2 + 32a + ^1 + AR^muj"^ - , 
7 = 7 = 1 + iV("-2)2 + 32a, 

where s = ^J{n — 1)2 + ^mEB? + AmR^uj'^ + 16a — 8&. The interval (0, 1) 9 t corresponds 
to the interval (0, 1) 3 z, therefore the requirement on asymptotic of the function f{t) near 
the point t = implies _ 

w{z) = F{a,l3;^]z). 

Also due to 

7-5-/3 = -iVl + 4i?4TOcj2 < 0, Re5>0 
and ((C.9(l . the requirement on asymptotic of the function f{t) near the point t = 1 implies 

P = ~k, fc = 0,l,2,... 

This leads to energy levels 

^'^^ 8^2 (^(4fc + 2 + - 2)2 + 32a) ^ - {n - i f ~ 16a + 8b+l 

+ ^ (4 + 2 + yF32FT32^) yiT^, . = , 1 , 2 , . . . 

Again multiplicities of these energy levels coincide with multiplicities of eigenvalues in 
propositions O El a-nd |5| 

The expression for radial eigenfunctions (up to an arbitrary constant nonzero factor) 

is 

rP+{r^-iy+ A i-iy (5)j- 4-'r2j 
(r.2 + l)^Pr+pi^'+iPr' ^0 J^k - jV- (7), (r2 + 1)2.- ' 

where p'"^\ P+'i P+°\ and 7 are given by above formulas for E = Ek- 
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7 Conclusion 

The possibility to find in explicit way some (but not all) eigenvalues for a Schrodinger op- 
erator characterizes so called quasi exactly solvable models In the present paper 
we have shown that the two-body problem on spheres S" with Coulomb and oscillator 
potentials is quasi exactly solvable for any n. A possible generalization for other compact 
two-point homogeneous spaces is an open problem. 

The quasi exactly solvability here is an attribute not of a radial differential equation 
(|32|l . but of the whole problem. It stems from the two causes. The first cause follows 
from the fact that we restrict our consideration on the subspace of C{G,K,iJ,) (see (|12|l ') 
consisting of common cigcnfunctions for operators Dq, Di, D2 and optionally D3. For 
every such eigenfunction one gets a separate radial differential equation (|32(l . For the 
Coulomb and oscillator potentials this equation can be reduced to the Heun one, but the 
further reduction to the hypergeometric equation using Maier's scheme is possible only for 
some eigenfunctions (just for those that satisfies the equation 1001)) and this is the second 
cause. 



A Orthogonal complex Lie algebras and their repre- 
sentations 

A.l Lie algebra !Bfc 

Here is a brief description of the simple complex Lie algebra 03^ = so(2fc + 1, C) (see 02], 
gH and gHI for details). 
Denote 








1 
1 



e GL(i), i G N. 



\1 ... 0/ 

Consider the Lie algebra Sfc = so{2k + 1, C) as 

= (Aegl(2A: + l,C)| S2k+i + S2k+iA ^ O) 



(A.l) 



Following I86j we shall enumerate the rows and columns of A g *8fc by the indices 
— fc, . . . , —1, 0, 1, . . . , /c. The convenience of such notations is due to the fact that sub- 
algebras Si C ?Bfe, i < k correspond to indices of rows and columns from —i to i. 
It can be easily shown that a matrix 



A^Y^ a^jEij G 0[(2fc + 1, C) 



belongs to *8fc iff 
secondary diagonal. 



0, which means that A is skew-symmetric w.r.t. its 



Let Fij — Eij — E^j,^i. It is easily seen that 



S F 



The algebra Q5fc is spanned by elements Fij with i > — j. Evidently, Fi^^i = and 

Elements Fa, i = l,...,fc form a base of the Cartan subalgebra t)k C Sfe, which 
consists of elements of the form 



X = dia.g{-Xk,-Xk-i, ■ ■ .,-xi,0,xi, . . .,Xk-i,Xk) ■ 
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Let Si e i)l such that ei{X) = Xi, i.e. is a base in dual to Fi^i, i — 1, . . . ,k. Define a 
symmetric nondegenerate biUnear form (•, •) on 25fe as 

(A,B)^^trAB, (A.2) 

which is proportional to the Killing form. Clearly, 

{Fij,Fgp) = Sip6jg, i > -j, q > -p. 

In particular, Fa, z = 1, . . . , fc is an orthogonal base in t)k- 

The form generates the isomorphism >€ : f)/c — > by the formula >c{X) = 

{X, •). Specifically, >c{Fi^i) = Si and Si, i = 1, . . . , fc is an orthonormal base in f)^. w.r.t. 
the form 

(/i,/2)* := {K-'ih),K-\f2)), fuh e ^)^ 

Using this notation one can describe the standard form of the root system for 5Bfc in 
the following way. Let 

$ajfc := (±£,,±e, ±ej| i ^ j, i,j = 1, . . . ,fc) 
be a root system in Q5fc, 

■= {£t,£'i + - Sjl i> 3, «,i = 1, ■ • ■ 
be a system of positive roots, and 

Asfc := (ai = £i, ai = £i- £-i-i| i = 2, . . . , fc) 

be a system of simple roots, corresponding to the inverse lexicographic order. A subalgebra 
QSj C Q3/c, i < k corresponds to root systems and A^b^- 

Let La be a root subspace in 05^, corresponding to a root a € •I'Sfc- Then 



L_e^ = span(i^oi), -^s; = span(Fio), ^£,-6^ = span(i^jj), 
^£,+6, = span(i^i = span(i^_ij), i,j = 1, . . . ,fc. 



(A.3) 



Fundamental weights for *Bfe are 



Ai = ^X^'^J' '^^^ =X!'^J' * = 2,...,fc. 

Let 

A = ^ A^Aj, A^' e Z+ := (0) U N 

be a dominant weight and V^(A) be an irreducible finite dimensional Q3fc-module with the 
highest weight A. All finite dimensional irreducible representations of 58^ are of this form, 
modules V{X) with different A are not isomorphic to each other, and V{X) corresponds 
to a (single valued) representation of the group S0(2fc + 1) iff Ai is even. The dominant 
weight A can be written in the form 

k 

A = ^ niiSi, nik ^ mfe_i > . . . > mi > 0, (A. 4) 
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where either all rrii £ Z+ or all rui G Z+ + i. Even values of Ai corresponds to nii E 
Let S be the sum of fundamental weights. Then it holds 



2 ^ V 2 

The universal Casimir operator C G ?7(*Bfe) is 
fc 

^ = E(^'^ + {^«'^o»})+ E + (A.6) 

1=1 i>i>0 

The following formulas are valid for any semisimple Lie algebra: 

C\y = {{6 + \,S + X)-{S,S)) id, (A.7) 

dimV{X) = ll{X + S,a) / ll{S,a), ( A.8) 

Q^O / ayO 

where a >- means a positive root. 

For any semisimple Lie algebra g and its Cartan subalgebra f) the module V{X) can be 
decomposed into the finite direct sum of weight subspaces 

^(A)=0F^(A), /iGf)*, 

where Vw £ V^(A),Vft, G f) it holds h{v) — fi{h)v and the sum is over weights of the form 

A - ^ iaa, ia e Z+. 

Besides, for any root a of g one has 

: V^iX) ^ V^+^iX), U£L^. (A.9) 



A. 2 Lie algebra T>k 

The Lie algebra Dk is the subalgebra of 55^, consisting of matrices whose column and rows 
with the index vanish. We shall discard these null row and column and shall enumerate 
other rows and columns oi A £ Dk by the indices — fc,...,— as before. The 
Cartan subalgebra \)k C Sfe is the same as in the ^Bfc-case. Describe the S)fc-case briefly, 
emphasizing differences from the Q3fc-case. 
Now one has 

:= {±e^±£j\ i ^ j, i,j 1,. . . ,fc) , 
<I>+^ := {£^+£■i,e^ -Ej] i > j , i , j = 1 , . . . , k) , 
^Vk := (ai = £i + £2,ai = £i - £i-i\ i = 2, . . . , fc) . 

The root subspaces L±i,.±^^ are the same as in *Bfe-case. 
Fundamental weights are 

Ai = ^E'^J' ^2 = + ^E'^J' = E^J' * = ^'---'^- 
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The sum of fundamental weights is 

k k 



2 

1=1 

A dominant weight 

fc 

A = ^ X^Xj, X^ e Z+ (0) U N 

now has the form 

fc 

A = ^TOiEi, mfc ^ rrik-i ^ . . . ^ m2 > |mi|, (A. 10) 

1=1 

where either mi E Z,mi E Z+, ?' ^ 2 or mi e Z + ^,771^ G Z+ + i, i ^ 2. Again 2)^- 
modules with integer nij, j ^ I, . . . , k correspond to (single valued) representations of the 
group S0(2fc). 

The universal Casimir operator C E U{'Z>k) is 

fe 

C^Y.^l + E m,,F,,} + {F,^^,,F^,,}). (A.ll) 

i=l j>j>0 

A. 3 Restrictions of and 2)fc-representations. 

The following results were found in 0^ (see also \bO\). 



Let V<8fc (A) be a simple *8fc-module with a highest weight (jA.4|l and V^^ (A) be a simple 
£)fe-niodule with a highest weight 

fc 

A' = ^m^e,, m'f^ ^ m'^.^ ^ . . . ^ > |m;|. 

i=l 

Proposition A.l. T/ie restriction V2jj.(A)|jj^ of the irreducible ''^k-representation onto 
any suhalgebra T)k C expands as follows 

^■B.WIs. -0^s.(A'), 

A' 

where the summation is over all A' such that 

nik ^ m'f. ^ ruk-i ^ . . . ^ m'2 ^ toi ^ m'^ ^ —mi 
and a/Z are integer or half integer simultaneously with mj . 
Let V<Sk-i{^') be a simple 5Bfc_i-module with a highest weight 
fc-i 

^' = E "^fe-l ^ ™fe-2 55 • • • > "4 m-i tfS 0. 

i=l 

Proposition A. 2. The restriction Vj)i^{X)\tg^ ^ of the irreducible Dk-representation onto 
any subalgebra 53fc_i C Sfc expands as follows 

A' 

where the summation is over all X' such that 

m-k > 'TiJ.-i ^ TOfc-i ^ . . . ^ m2 ^ JTi'j > |mi| 
and all m'^ are integer or half integer simultaneously with mj . 
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B Self-adjointness of Schrodinger operators on Rie- 
mannian spaces 

Here we shall formulate two results concerning the self-adjointness of Schrodinger operators 
on Riemannian spaces, which is used in sectional 

The first theorem is a result from restricted onto the scalar case. 

Theorem B.l. Let M be a Riemannian manifold of a bounded geometry, dimM = £, 
and II be the measure on M generated by its metric. Suppose also that the potential V can 
be represented in the form V ~ Vi +V2, where real valued functions Vi, V2 are as follows: 
Os^Vi e Cl^{M, fi), ^ ^2 e CP{M, /i) for p = £/2 if £ ^ 3, for p > 1 if £ = 2, and for 
p^l if£^l. 

Then the operator Hy = — A +V is self-adjoint with the domain: 

Dom{Hv) ^ (^ue W^^^{M,fi)\ J Vi\u\^dfi < +00, Hyu G ^^(m,^)^ , (B.l) 

where Hyu is understood in the sense of distributions. Here W^''^{M, fi) is the Sobolev 
space, consisting of functions on M that are in >C^(M, /i) with their first derivatives. 
Also Vu e Cl^{M,ii) for u e Bom{Hv). 

The definition of a Riemannian manifold of a bounded geometry can be found in 
[52|. Note that compact and homogeneous Riemannian manifold is always of a bounded 
geometry. 

If the potential V is not in Cl^^{M",iJ,) then theorem IB. II is not applicable. If instead 
V is bounded from below, one can try to restrict the Schrodinger operator onto some 
submanifold M' of such that V\J^.J, € Cl^^{M' , fi) and construct the Friedrichs self- 
adjoin extension of — A +V from the initial domain C^{M'). This procedure is 
physically motivated for instance in the case when V —^ +00 near the boundary of M' 
and therefore wave functions should vanish near this boundary. 

Let us turn to the accurate mathematical description. Let M' be an open connected 
submanifold of a Riemannian space of dimension £ with a metric g and an induced 
measure fi. We do not suppose that M' is complete w.r.t. the Riemannian structure 
induced by the Riemannian structure on M^. Let y ^ C G M be a real valued function 
from Cl^^{M' , fi) and H' ~ — A +V be a Schrodinger operator with the domain C^(M'), 
consisting of all infinitely smooth complex valued functions in M' with compact supports. 
Not loosing generality we suppose that C — 1. Let Hp ^ id be the abstract Friedrichs 
extension of H' [^Sl- We need a precise description of 'Dom{Hp). 

The operator H' generates sesquilinear nonnegative form qh' by the equality 

QH' {(p, ^) = {H'lp) ^dfi 

JM' 

with the domain C^{M'). Evidently, its closure is 



JM' 

with Bom{qHp) C C'^{M',fi) being a closure of C^(M') w.r.t. the inner product HB.2|I . 
where V is the gradient operator given in local coordinates by the equality 

The operator Hp is defined by the identity 

/ {g{W<f,'Vi(;) + V<fiip)dfi = / ipHpiljdn, Vip e I)om{qHp),ilj E BomiHp). 

JM' JM' 
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Thus 

Hf^P = a + V^b)^,^, , ^ e BomiRp). (B.3) 
Theorem B.2. The domain of the operator Hp is 

(V e W'^' (M',^) I e Aoc(m',m); (- A ^ + FV)dist e -C^C^^'^m)) 



and Hp acts by formula (|B 

The proof of this theorem repeats mutatis mutandis the proof of theorem X.27 from 
[SS! using the generahzation of the Kato inequality for Riemannian spaces jSH- 

C Some Fuchsian differential equations 

For convenience of references we collected here basic facts concerning some Fuchsian dif- 
ferential equations: the Riemannian equation and the reducibility of the Heun equation 
to the hypergeometric one. 

The linear differential equation 

w(")(z) +pi(z)w("-l)(z) + . . .+pn{z)w{z) = (C.l) 

on the Riemannian sphere C = P^(C) with meromorphic coefhcients Pi{z), i — 1, . . . ,n is 
Fuchsian if for any € C its solutions has no more than a power growth as z tends to 
zq in some cone^, not containing a neighborhood of its vertex zq. A point zq is regular for 
this differential equation if all solutions of (|C.1|I are holomorphic in some neighborhood of 
Zq; otherwise zq is a singular point. 

It is known 55 that equation IC.ip is Fuchsian one iff 

PH^) = — 



n i^^-^kY 

k=l 



for some finite potentially singular points zi, . . . , Zm € C and polynomials qi{z) of degrees 
i{m — 1). One can find characteristic exponents p^^*' of (|C.ip at the point Zk by 
the substitution w{z) — (z — z^Y'' into 1C.1(I and keeping only leading terms as z — > 
Zfe. This procedure gives an algebraic equation of the n-th degree for p^^''\ Denote by 
^(2fc)^ I = 1, . . . ,n its solutions for all points Zk, k — 1, . . . ,to. The substitution w{z) — 
z~P^ ' similarly gives characteristic exponents p^^^ , ■ • • , Pn^^ in the infinity. Characteristic 
exponents satisfy the Fuchs identity: 

n m+1 -J 

k=i 

where p] ■= pI ■ 

One can find characteristic exponents also for a regular point. If a point z is regular, 
then characteristic exponents for this point are 0,l,...,n— 1. The sufficient condition for 
the regularity of z is the regularity of coefficients Pi{z), i = 1, . . . , n at this point. 

An information on singular points and corresponding characteristic exponents of equa- 
tion can be encoded in the Riemann P-symbol P{A;z}, where the first row of a 
matrix A consists of singular points and other rows of A consist of corresponding charac- 
teristic exponents. 

®In a neighborhood of the infinite point one should use the local coordinate ( = 1/z instead oi z — zq. 
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Equation ljC.l|) of the second order with three singular points is called the Riemannian 
equation. Coefficients of the Riemann equation are completely defined by its characteristic 
exponents. Equivalently, the Riemann equation is completely defined by its P-symbol. In 
this case the Fuchs identity looks like 

fc=i 

and there are only five independent characteristic values. 

There are two types of variable change, transforming any Fuchsian equation into an- 
other Fuchsian equation. The first one is a linear-fractional (Mobius) transformation of 
the independent variable: 

z^t, z = ^^^,aS-Pj^O. (C.2) 

By such transformation one can move three singular points into three arbitrary points of 
C with the same characteristic exponents. 

The second one is a linear transformation of the dependent variable 

fz-ziY 

w{z) wi{z) = I w(z), (C.3) 

\z-z2J 

which conserves singular points, but changes the characteristic exponents 

(zi) (zi) (Z2) (22) ■ 1 n 

Using these transformation for the Riemannian equation one can move three singular 
points into the triple (0, 1, 00) such that pf^ = p'p = 0. If one denote p'f^'' — a, P2°°^ — (3 
and p^2^ = 1 — 7, then the Fuchs identity for this equation gives p^2^ — ^ ^ a — (3 that 
corresponds to the hypergeometric or Gauss equation: 

z(l - z)w"{z) + (7 - (a + /3 + l)z)w'{z) - al3w{z) = 0. (C.4) 

The P-symbol of equation ljC.4(l is 

1 00 

P{ Q a ; 

1 — 7 7 — a — /? (i 

Many quantum mechanical problems for constant curvature spaces can be reduced to 
this equation, while their Euclidean counterparts lead to its limiting cases, obtained from 
()C.4ll by confluence of singular points (such equations are not Fuchsian). 

We shall consider only solutions of (|(].4|l in the case 7 7^ — m, m e N. Solutions of 
(jC.lll . corresponding to different characteristic exponents near some singular point are 
called canonical solutions near that point. The series 



F(a,/3;7;z):=^^^)^^^, |z| < 1 (C.5) 

n=0 



where (a)„ := a(a -I- 1) . . . (a + n — 1), (a)o '■— 1, is the canonical solution of (|C.4|I . cor- 
responding to the characteristic exponent p^^^ = 0. The function F(a,(3;^] z), defined by 
(l(T5ll for \z\ < I can be analytically continued for z G C\(l, +00) [21], |57| . 

Evidently F{a, (3; 7; z) = F{P, a; 7; z). If a = —to 01 f3 = —m, to = 0, 1, 2, . . ., then 
F{a, (3; 7; z) is a polynomial of degree to. 
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Another canonical solution of ljC.4|) . corresponding to the characteristic exponent 
p^2^ = 1 - 7 for 7 ^ N, is 

z^'^F{a - 7 + 1, /3 - 7 + 1; 2 - 7; z). 
Canonical solutions near the singular point 2 = 1 are 

F{a,l3;a + (3--/ + z) 

and if7 — a — /3^Z also 

(1 - z)''-"-^i^(7 - a, 7 - /?; 7 - a - /3 + 1; 1 - z). 
Near the singular point z — 00 canonical solutions are 

z-'^F{a, a - 7 + 1; a - /? + 1; i), z-f^FiP, /3 ~ 7 + 1; /3 - a + 1; i) 

if a — /? ^ Z. If a — /9 G Z, then only one of these expressions is a canonical solution: the 
first if a — /3 > or the second if a — /3 < 0. 

There are expansions of F(a, z) through canonical solutions near the singular 
points z — 1 and z — 00 |57j . j58| . important for spectral problems. The first one, used in 
the present paper, is 

r(y)T(y -a - 3) 

F{a, /3; 7; z) = "-/^ ^\^. ^F(a, p;a + (3 - ^ + 1;1 - z) (C.6) 

r(7-a)r(7-p) 

+ ^^^w",X!!r^^ (^ - z)^— ^F(7 - a, 7 - /3, 7 - « - /3 + 1, 1 - I arg(l - z)| < ^ 

r(a)r(/3) 

if 7 — a — /? ^ Z. Here F is the gamma-function. It has no zeros and has poles of the first 
order at the points z = —m, m = 0, 1, 2, . . .. Its logarithmic derivative V'r(-z) := r'(z)/r(z) 
also has poles of the first order at the same points. 

For 7 — a — /3 e Z every summand at the right hand side of (|(].6|l is singular and it 
holds for TO = 0, 1, 2, . . . 

" + + r(a + TO)r(/3 + to) ^ n!(l - to)„ ^ ^ 

1 (all p '^-^ n'An + my. 

X (ln(l — z) — ^r(^ + 1) ~ ^r('T' + to + 1) + tjjrio^ + n + to) + ipriP + n + to)), 

(_l)rnr(a + /3-TO) ^ («)„(/?)„ 

F(a - to)F(/3 - to) ^n!(n + TO)!*^ ^ ' 

X (ln(l - z) - ipr{n + 1) - + to + 1) + ^ria + n) + ?Ar(/3 + «)), 
|arg(l - z)| < TT, |1 - z| < 1. 

In the case Re(7 — a — /?) < formulas (|C.6(I - (|C.8|I imply 

lim^F(a,/3;7;.)(l- z)-'^- = "^^y,^" ^^^^^ (C.9) 
z^i r(a)F(/3) 
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The Fuchsian equation l|C.l|) with four singular points by transformations ljC.2|) and 
(|C.3p can be reduced to the Heun equation 

where 0, 1, d, oo are its four singular points {d ^ 0, 1, oo) and a + /3 — 7 — (5 — e + l = 0. 
The corresponding P-symbol is 

To 1 d c» 
P< a ;t 
[1-7 l-S 1-e P 

Note that the accessory parameter q does not arise in this P-symbol. 

The theory of the Heun equation is much less explicit than the theory of the Rie- 
mannian equation. In particular, there are no explicit expressions of canonical solutions 
near different singular points through each other. Therefore there are only approximate 
methods for solving spectral problems connected with the Heun equation, using continued 
fractions (see for example [^2] and references therein). 

The substitution z = P{t) for a rational function P transforms equation IjCljl into 
another Fuchsian equation with generally a greater number of singular points. Therefore 
sometimes the inverse transformation can decrease the number of singular points of a 
Fuchsian equation.^ 

At the present time there is no a general theory of such reduction. However in 
there were classified all Heun equations IjClOjl that can be obtained by a substitution 
z = P(t) from the hypergeometric one ljC.4|) . By the inverse transformation these Heun 
equations are reduced to hypergeometric equations. 

The first condition for existing such reduction is the position of the point d. Let 

{Zl - Z^){Z2 - Z4) 



izi,Z2,Z3,Zi)c, 



{Zl - Zi){z2 - Z3) 



be the cross-ratio of four pairwise distinct points from C. It is well known that a 
cross-ration is invariant under Mobius transformations. The group S4, permuting points 
Zl, 22,2:3 and 24, acts on their cross-ration. The cross-ration orbit (^©^(s) ofs := (21,22,23, 
Z4)c.r. consists of points s, 1 — s, 1/s, 1/(1 — s), s/(s — 1), (s — l)/s € C. 

In general position this orbit consists of six points, but there are two exceptional 
cases: the orbit — l,i,2 and the orbit ^ ± -^i. If (21, 22, 23, 24)c.r. G then 

(21,22,23,24) is a harmonic quadruple. If (21, 22, 23, 24)c.r. — \ ^ "^ij then (21,22,23,24) 
is an equianharmonic quadruple. 

Points of a harmonic quadruple lie on a circle or on a line. By a Mobius transformation 
they can be mapped into vertices of a square in C. If (21, 22, ^3, 00) is a harmonic quadruple, 
then (21, 22, 23) are collinear, equally spaced points. If (21, 22, 23, 00) is an equianharmonic 
quadruple, then (21, Z2, 23) are vertices of an equilateral triangle in C. 

Theorem C.l ([60 ). All cases, when nontrivial Heun equation (|C.10|I (i.e. a/3 ^ or 
q 0) can be obtained from the hypergeometric one (|C.4|I by the rational substitution 
2 = P{t), are as follows. 

1. Harmonic case: d £ 0(s^{2). Suppose d — 2,^ then q/{a(3) must be equal 1, and 
characteristic exponents of points t = and t = d — 2 must be the same, i.e. 7 — e. 
The function P{t) is a degree-2 polynomial and can be chosen as P(t) = t{2 — t) — 
1 - (t - 1)^ It maps t = Q,2 to z = Q and t ^ I to z ^ l.^ 

^Generally, the inverse transformation does not conserve the Fuchs class of differential equations. 
*If d £ Oq^{s), then the quadruple (0, 00) can be mapped into the quadruple (0, l,s,oo) by a 
Mobius transformation, which transforms also parameters of equation IC.lOi . 
^This transformation was found already in Itill . 
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// additionally 1 — 6 = 2(1 — 7), then P{t) can be chosen also as degree-A polynomial 
4(t(2 -t)- which maps t = 0, 1, 2 to z = 1. 

2. d €z 064(4). Suppose d = A, then q/{aP) must be equal 1, characteristic exponents 
of points t = 1 must be double those of the point t = d = A, i.e. 1 — 6 = 2(1 — e), and 
t = must have characteristic exponents 0, 1/2, i.e. 7 = i. The function P{t) is a 
degree-?) polynomial and can be chosen as [t — 1)^(1 — j). It maps t = to z = 1 
and t = 1,A to z = 0. 

3. Equianharmonic case: d £ ©64(5 + "^i)- Characteristic exponents of points t = 
0,1, d are the same, i.e. j = 6 = e. Suppose d = ^ + -^i, then q/{af3) must be 
equal ^ + -^i- The function P{t) is a degree-? polynomial and can he chosen as 

[1 — t/(^ + -^i)^ • It maps t = Q,l,d to z = 1 and t = q/{a[3) to z = 0, thus 
creating a new singular point. 

If additionally j = 6 = e = |, then P{t) can be chosen also as degree-Q polynomial 

which maps t = 0, l,d,q/{af3) to z = 1. 

4. d <E OeAh + ^i)- Suppose d = ^ + ^i, then ql{a(3) must be equal \ + ^i, 
characteristic exponents of the point t = d must be 0, 1/3, i.e. e = 2/3, and points 
t = 0,1 must have characteristic exponents 0,1/2, i.e. ^ = 6 = 1/2. The function 
P{t) is a degree-A polynomial and can be chosen as 




It maps t = Q.l to z = 1 and t = d,q/ {af3) to z = 0. 

5. de Oe^{\ + ^^^i). Suppose d=\ + "g^ i, then q/{aP) must be equal \ + 

characteristic exponents of the point t = d must be 0, 1/2, i.e. e = 1/2, and points 
< = 0, 1 must have characteristic exponents 0,1/3, i.e. 7 = 5 = 2/3. The function 
P{t) is a degree-5 polynomial and can be chosen as 

.2025\/T5 , , / 1 VlE.y 

-i — t{t-l)\t — i . 

64 ^ \ 2 / 

It maps t = 0,1, q/{aP) to z = and t = d to z = 1. 

Note that there are three independent parameters in the first case of theorem (jC.l|) (for 
example: a, /3, 7) and all other cases contain only one or two free parameters. It means 
that the first case is more rife in applications. In fact, it is the only one, which occurs in 
the present paper. 

D The proof of two expansions 



Expansions H25(l . H31(l were obtained in PH] and "^S* using the theory of Yangians, which 
is a part of the quantum algebra. Here we give an independent proof of these expansions 
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from a classical point of view using one result from . The initial idea of this proof is 
due to A.I. Molev. 

Let A be an associative algebra over C, generated by elements Z+, Z_, F and relations 
[F, Z+] = 2Z+, [F, Z_] = -2Z_, [Z+, Z_] = -\f^ + qF, q G C. (D.l) 

Let also t : A ^ Homc(V^) be its irreducible linear representation in a finite dimensional 
complex linear space V . An arbitrary linear operator in V has at least one eigenvalue. Let 
Fx = r]x, rj & C, X ^ V^X 0, then relations (|D.1|I imply FZ+x = iv + 2)2'+x, FZ^x = 
(77 — 2)Z_x. Since dime V < 00 there is a vector £ V, ^ such that Fvi, = 

ly, Z+Vi, = 0, e C. Let iV-2j := Z-Vi,, j G Z+, then 

Fu^ = 771;^, Z^v,j = v,j^2, V?/ G Li, := - 2Z+. (D.2) 
Let /Lt be a root of equation 

fi^ + + + 2iy - Aq = 0. (D.3) 

Lemma 3. It holds 

Z+Vr, - ^(?7 - /x)(r? - i/)(r? + /i + 2)(?7 + 1/ + 2)v^+2, V & L^- (D-4) 



Proof. For r] — v equality l|D.4l) is obvious. Let a € L^, a ^ v and suppose that l|D.4|l is 
valid for any rj €z such that rj > a. Then one gets 

Z+Va = Z+Z^Va+2 = [Z+, ZJ\Va+2 + Z_Z+Va+2 = (^^^^'^ + 

+ -^rCa + 2 - + 2 - + M + 4)(a + + A)Z^Va+i 
16 

+ 2)^ + q{a + 2) + -^(a + 2 - ^)(a + 2 - + + 4)(a + v + A)\vo, 
z lb 

-^(a - /i)(Q; - i/)(q; + /i + 2){a + v + 2)va+2, 
lb 



+2 



due to (|D.1|I . (|D.2(I . I|D.4|I and the following identity 

(a + 2 - ^) (a + 2 - I/) (a + /i + 4) (a + J/ + 4) - (a - /i) (a - (a + ^ + 2) (a + i/ + 2) 

= 8(a + 2)3 - 16g(a + 2). 

This completes the induction. □ 



Due to (jD.2|) nonzero vectors u^, e are linear independent, therefore from dime ^ < 
oo one gets ^ 0, = for some vi £ L,^. Since the representation t is irreducible 

one gets 

V = span{v^^,v^^+2, ■ • ■ , dime = -{ly - i^i) + 1 (D.5) 

due to HD.2|I and (jD.4|) . Equation HD.4|I for rj = vi - 2 implies 

(j/i - 2-^)(i^i = 0. (D.6) 

Let now Z+ = 13+ , Z_ = 13^, F ~ Fkk, k ^ 2 ior operators D^, , Fkk from section 
15.11 V :— V<Skimk£k + "^fe-iE/c-i) and therefore 

q=l {ml + ml_^ + {2k - l)mk + {2k - 3)mfc_i) + ^{2k - l)(2fc - 3). (D.7) 
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Due to theorem 9.1.12 from 62 the ^-module V is irreducible, therefore lemma |21 is 
applicable. We shall demonstrate that 

(i/i-2-/i)(z/i+/i)^0. (D.8) 

Indeed the y^-module V^s^irnkSk + "mk-iSk-i) does not contain weights jsk for \j\ > ruk 
P^ . therefore 

(i^i-l)2 + (,y + l)2 ^2(mfe + l)2. (D.9) 
If (z^i 2 — = then — 1| = |/i + l| and equation l|D.3|l leads to the inequality 
2 + iq={fi + lf + {i^ + 1)2 <C 2(mfc + 1)2, 

which is equivalent to 

m^._i + {2k - 3)(mfc + nik-i) + ^{2k - l)(2fc - 3) (D.IO) 

due to HD.7|I . Obviously inequality (|D.10(I is impossible. 

Thus from (|D.6|) . I|D.8(I and (|19|l one gets z/ = —vi = nik — nik-i- Now expansion H25(l 
follows from l|D.2(l and l)D.5|l . 

The consideration for the space := Vj)j.(mfcefe + mfe_iefc_i), fc ^ 2 from section 1^?^ 
is similar. Here 

<? = ^ ("ifc + mg_i + 2{k - l)mfe + 2{k - 2)mfe_i) + (fc - l)(fc - 2) 

and the conjecture (z^i — 2 — /i)(i^i + /i) = now implies due to (|D.9|I 

2(A: - 2)mk + (mfe_i + fc - 2)^ + A:(A: - 2) 

that leads to fc = 2, nik-i — and to the equality in l)D.9(l . This yields 

V = -vi = nik - \mk.^i\. (D.ll) 

The last possibility I'l = -ly in (|D.6(I also leads to l|D.ll() due to and (|D.5(I . Now 
expansion (|31|l is a consequence of ljD.2p . (|D.5|) and (jD.ll|l . 
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